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Abstract: We investigate charged dilatonic black p-branes smeared on a transverse cir- 
cle. The system can be reduced to neutral vacuum black branes, and we perform static 
perturbations for the reduced system to construct non-uniform solutions. At each order a 
single master equation is derived, and the Gregory-Laflamme critical wavelength is deter- 
mined. Based on the non-uniform solutions, we discuss thermodynamic properties of this 
system and argue that in a microcanonical ensemble the non-uniform smeared branes are 
entropically disfavored even near the extremality, if the spacetime dimension is D < 13 +p, 
which is the critical dimension for the vacuum case. However, the critical dimension is not 
universal. In a canonical ensemble the vacuum non-uniform black branes are thermody- 
namically favorable at D > 12 +p, whereas the non-uniform smeared branes are favorable 
at D > 14 + p near the extremality. 
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1. Introduction 

Resolving the structure of vacuum solutions of Kaluza-Klein theory and finding the final 
fate of the Gregory-Laflamme (GL) instability |l|, ||] of black objects has been an intriguing 
subject in the last decade. Recent progress in the subject has revealed many aspects of the 
problems, but at the same time several important issues remain to be resolved ||, 0]. 

Perturbative analysis of the vacuum system on a circle has been made extensively, 
in particular for black holes |5], ^, [?], || recently. For uniform black strings, perturbative 
non-uniform deformation is allowed only for the critical uniform string that has the static 
GL mode ||. However, the most interesting state of black objects on a circle is in the 
regime beyond perturbations, and so far only a numerical approach has succeeded. The 
fully nonlinear static solutions with nontrivial horizon geometry have been constructed in 
6 dimensions (and, partially, 5 dimensions) Jl^, |ll], |l2|, 13, 14]. The numerical solutions 



are constructed by employing the conformal gauge that is a general metric ansatz under 
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axisymmetry and advantageous to numerics. The constructed two branches of a vacuum 



state, black hole and black string, help us to make a phase diagram in this theory [14]. 

To push forward our understanding of vacuum states in many directions, it might be 
useful to employ a metric ansatz proposed by Harmark and Obers (HO) |1EJ, which involves 
just two undetermined functions. Although the metric ansatz was originally hypothetical 
because its proposers could not provide a proof of self-consistency, it has been shown 
that their simple ansatz is equivalent to the general conformal ansatz by an appropriate 



coordinate transformation |16|, 17, 18]. A key point is that one can also show the consistency 
of the associated boundary conditions. 

While vacuum solutions represent a rich vacuum structure, charged black objects on 
a circle are also an interesting subject. A system of charged black objects has in general a 
thermodynamically stable parameter regime, and related to this fact, there is a conjecture 
that correlates dynamical stability and thermodynamic stability. The conjecture is known 



as the Gubser-Mitra (GM) conjecture [|T9|, ^(J (or the correlated stability conjecture j2l|]), 



which states that for systems with a translational symmetry and an infinite extent the 
dynamical GL instabilities arise precisely when the system is thermodynamically unstable. 
There is a lot of evidence for this conjecture, and no fully demonstrated counterexamples 
have been discovered so far (see, e.g., [21, 22, |2^]). Another aspect of introducing charges 



into the system on a circle is connections with non-gravitational theories [p4| , 25]. Besides, 
existence of charged non-uniform black branes having higher entropy than their uniform 
counterparts is conjectured in Ref. p6fl . 

In this paper, we investigate charged dilatonic black branes smeared on a transverse 
circle, which we call smeared black branes. Following the perturbation theory developed 
by Gubser Q , we study the stability and possible non- uniform deformation of the smeared 
black branes. To perform the perturbations, we employ the HO metric ansatz, which 
greatly simplifies the analysis in many respects. Based on the solutions of non-uniform 
smeared black branes, which include vacuum case, we study the thermodynamic properties 
of the system in detail. We think that the perturbative analysis developed in this paper will 
give a good theoretical basis for constructing fully nonlinear solutions of the HO ansatz, 
as the original perturbative solution played an important role in constructing nonlinear 
solutions. 

The organization of this paper is as follows. In the next section, we review the smeared 
black branes in the HO metric ansatz and derive thermodynamical quantities. In Sec. |||, 
we perform a perturbative analysis up to the third order. The thermodynamics of the 
smeared black branes is discussed in Sec. |], and Sec. [| is devoted to a conclusion. 

2. Smeared black branes 
2.1 Harmark-Obers ansatz 

The part of the classical supergravity action relevant for the considerations in this paper is 



Id = — %r~ / cPxy/^g R - -(<90) 2 - , 1 N , e^F 2 2 



(2.1) 
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where F p+ 2 is a (p + 2)-form field and is a dilaton field. The coupling constant a is 
a 2 = 4 — 2(p + l)(d — 2)/(D — 2) with D = p + d + 1. From this action, we obtain the 
following equations: 



□oi - a r a<t> F 2 
^~2(p + 2)\ e ^+ 2 ' 



(P + 2)F„ 



F 2 



(2.2) 



a. 



-ge a ^F^^ 2 '"^ p+2 ) = 



The form field must satisfy the Bianchi identity, d^F^ ...^ +2 ] = 0. 

The ansatz proposed in |15| ] for a charged dilatonic black p-brane with transverse space 
R 6 ^ 1 x S 1 is 



ds 2 D 



d-2 



-fdt 2 + £ (dx*) 2 + ^ (^) 2 (y^ 2 + ^ 2 + 5i? 2 d0 2 d _ 2 



/ = l- 



d— 3 



e 2 * = H a , 



H = l + 



Rq 

~R 



d-3 



sinh 2 a, 



(1 - fl-^cotha, 



(2.3) 



where A = A(R, v) and B = B(R, v) are undetermined functions, and A is the potential of 
the form field F = dA. The constant a is a parameter of electric charge. The coordinate 
v has the periodicity v = v + 2tt, and the asymptotic size of the circle is L. The spacetime 
has an event horizon at R = Rq. The uniform smeared black p-brane is given by setting 
A = B = 1. The consistency between this ansatz and the conformal form for the vacuum 
case and the charged non-dilatonic case are given in [16, 17] and ][[q| , respectively. The 
consistency for the most general (charged dilatonic) case is discussed in |J7) by applying a 
Lorentz boost and U-duality. 

The remarkable character of this metric is that the EOMs for the two unknown func- 
tions A and B are independent of the value of p and the charge parameter a. Thus the 
EOMs are the same as those of the neutral non-dilatonic black branes on R^ 1 x S , which 
we represent in Appendix |A]. Consequently, the problem of finding solutions of charged 
non-uniform smeared branes is mapped to the problem of neutral non-uniform strings with 
the effective spacetime dimension T>, 1 



D = D -p = d+ 1. 



(2.4) 



For example, the non-dilatonic charged cases are (D,p, d; a) = (11, 5, 5; 0) and (10, 3, 6; 0), 
and their effective dimensions are D = 6 and 7, respectively. The two cases correspond to 
M5- and D3-branes. In this paper, we take D as a parameter, and we study perturbations 
in a wide range of D. 



1 Generating (near-extremal) charged solutions from uncharged solutions via the M-theory lift-boost- 
reduction is possible. This kind of solution-generating technique was used in [E4l G5l |2q]. 
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2.2 Thermodynamic quantities 



Thermodynamic quantities of the smeared black branes are calculated from the metric 
(|2.3j). The mass, temperature, entropy and charge are given by 



M 
T 



( LRr 



d-3 



V 2vr 
d-3 



(d-3) 



d-2 
d-3 



X + sinh a 



2RqL\/ Ah cosh a 



S = iirSVAh 



Q 



( LRr 



V 2?r 



( LRq 

d-3 



d-2 



cosh a, 



(d — 3) cosh a sinh a, 



(2.5) 



where we have used S = LV^rirf_2/16vrGD- ^4h(i?o) = ^4(i?o ; ^) and x = x(-^o) are con- 
stants for a given horizon radius. The latter is defined by the asymptotic behavior of the 
metric functions: 



A,B~1- X 



~R 



d-3 



for > Rq. 



(2.6) 



^ is the volume of x l (i = 1, 2, • ■ • ,p) directions. It is not necessary to take the coordinate 
x l to be periodic, but we do it formally so that we may present a finite expression of mass 
and charge. The volume of the unit n-dimensional sphere is f2 n = 2ir( n+1 ^ 2 /T [(n + 2)/2]. 
For a regular event horizon, is independent of v so that the zeroth law of thermody- 
namics holds. The tensions, which are asymptotic charges accompanied by the S 1 com- 
pactification, are given by p7|, 29, ^0|, 31] 



LT 



V 2tt 



[l-(d-3)(d-2) X ] 



LiTi = 9 {^r) d 3 [ 1 + ( d - 3 )( sinh2 «-x)] 



(2.7) 



where T and Tj are the tensions along the u-direction and along the i-th world-volume with 
circle length Lj, respectively. Defining the dimensionless tension n by 



LT 



n 



M — Mgi 



(2.8) 



where M e i = [(d — l)LjTj — LT — M] /(d — 2) is "electric" mass [25|, the Smarr formula 
then takes the form 



TS 



d-2-n 
d-1 



(M — M e l). 



(2. 



Note that for the neutral case (a = 0), we have M e \ = and that then the dimensionless 
tension becomes the one commonly employed for vacuum black branes. 
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The specific heat Cq and the isothermal permittivity €t are computed as 
dM\ V p Sl d _ 2 L d - 1 (R \ d ~ 2 [(f-3 + (d-l)cosh2a]cosha 



° Q ~ \ 8T ) Q AG \2tt) -(d-3) + (d-5)cosh2a ' (2 " 10) 

' / = iS) = Vp ^:t d ~ 2 ( ^ (d - 3) l(d - 3) - (d - 5) cosh 2a] cosh 2 a, 



d$J T 32vrG \2-k J 



where 



$ = tanha (2.11) 



is an electric potential energy at the horizon. The specific heat changes the sign from 
minus to plus at a critical value of a, if > 6. For a fixed charge Q, this implies that 
a thermodynamically stable region exists for the smeared black branes in the parameter 
region in which black branes suffer from dynamical GL instability pq| . However, we should 
consider the thermodynamic ensemble that allows the charge to change. In other words, 
we should also take into account the isothermal permittivity, which probes the thermo- 
dynamic stability under changes of the charge |23|, [3^]. One finds that the specific heat 
and isothermal permittivity always have opposite signs from the above expression. Conse- 
quently, the two conditions Cq > and ex > 0, which guarantee thermodynamic stability 
in the grand canonical ensemble, cannot hold simultaneously. That is to say, the smeared 
black brane background is never thermodynamically stable. This is perfectly compatible 
with the existence of a static mode that is a signal of the onset of unstable modes. We 
discuss the static mode in detail in the next section. 2 

3. Static perturbation 

According to the perturbation theory in Ref. H, we perform static perturbations of the 
smeared black p-branes. 3 In general, the static perturbations are easier to perform than 
dynamical perturbations, and they allows us to construct a non-uniform solution as well 
as to determine the GL critical waves. For convenience, we transform the metric functions 
and the coordinates as 

A = e a , B = e b , 

x = —v, y = —R. (3.1) 
2vr 2vr y J 

It is possible to rescale the entire metric so that the event horizon can be located at y = 1 
in the new coordinates. In the following analysis we adopt this normalization. 

We expand the metric function X(x,y) (X = a,b) around the uniform solution as 
follows: 

oo 

X(x, y) = J2 e " x «(y) cos nKx - ( 3 - 2 ) 

n=0 



2 A partial preliminary analysis of dynamical s-wave perturbations is found in 

3 Readers who are interested in the details in this section can examine a Mathematica notebook 
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Here, e is a small parameter of expansion. The function X n and wave number K = 2n/L 
are also expanded as 

oo oo 

X n (y) = Y,e 2p X n , p (y), K = Y,e 2q k g , (3.3) 

p=0 9=0 

where Xoo(y) = is imposed. Some leading order terms are 

X(x,y) = eX lfi (y)cosk x + e 2 [X 0jl (y) + X 2 ,o(y) cos 2fc a;] 

+e 3 [Xx,i(y) cos /co^ — Xifi(y)k\x sin fcox + A^^y) cos Sk^x] + 0(e 4 ). (3.4) 

Having expanded the metric functions as above, all we have to do is to solve the 
ordinary differential equations (ODEs) for X„ jP (y) at each order of e. Before entering into 
detail, however, it would be suitable to outline our calculations. First of all we should 
notice that in this perturbation theory the asymptotic circle length (or the wave number) 
is not fixed at all, and we have to determine it by solving the perturbations. At the first 
order O(e), we look for a solution X^o, which corresponds to the GL static mode. Thus, we 
can show the instability of the uniform solution by finding the GL critical wave number k^. 
At the second order 0(e 2 ), a back-reaction leads to nontrivial X q i and ^2,0- O ne of them 
is the Kaluza-Klein (KK) massive mode, which falls off exponentially for a large y. The 
other is a massless mode, which determines the mass of smeared branes as an asymptotic 
charge. Therefore, one might consider the second-order calculation to be enough to see a 
thermodynamical symptom of the non-uniform deformation. It is known, however, that to 
calculate entropy difference between non-uniform and uniform solutions for the same mass 
one needs third-order perturbations, at least in the uncharged case. This is also the case 
for our general situation. 

It is important to note the advantage of working in the HO metric. The metric ansatz 
reduces the number of equations of motion (EOMs) greatly. As discussed in the original 
paper , the EOMs do not contain the charge parameter a and the spatial world- volume 
p. The EOMs for most general charged diatonic cases are equivalent to those for neutral 
black strings on a cylinder R^ 2 x S 1 with T> = D -p = d+1. That is, we do not need 
to perturb the form and dilation fields at any order. 

It is worthwhile to consider the meaning of the parameter e. In a leading order, this 
parameter e can be related to the parameter used in the literature Q, 

x l(R^ L _ l \ (35) 

Here, R ma x and R m in are the maximum and minimum of the proper horizon radius of 
the black object so that A means geometrically the non-uniformity of the black brane. 
Substituting the expansion ( |3.2| ) into the metric ( |2.3| ), the parameter e is identical to A at 
the linear order if we set &i,o(l) = 2, which is possible due to the linearity of the differential 
equation at this order. Consequently, the expansion parameter e describes in the above 
sense how the black brane is non-uniform. However, the relation between e and A becomes 
nonlinear when we take into account higher-order perturbations. Because of the nonlinear 



- 6 - 



corrections, the naive identification between e and A is apparently invalid in the regime. 
After all, e is just an expansion parameter and is no more and no less than the parameter 
controlling the amplitude of the perturbations. 

3.1 First-order perturbation Xi t a 

Let us begin with the first-order perturbation. Substituting the expansion ( |3.4|) into the 
Einstein equation (|2,2|), we obtain a master equation for 61,0, which determines the GL 
critical wave number ko, 



£.161,0 = 0, (3.6) 



where 



= (d-l)f-(d-3)(3d-8)f + (d-3) 2 
[V) yf[d-3 + (d-l)f] 

V ^ iy) = yifld-3 + (d-l)f] • (3 ' 7) 

After some manipulation, one recognizes that the problem is the same as solving the 
Schrodinger-type equation with energy — although the potential is rather complicated 
in the present case. 4 

One can solve algebraically for ai o, with the result 

2(rf-2)/(6 1 , + 2/6' 1|0 ) 
ai < o(y) = d-3 + (d-l)f ' (3 ' 8) 

where the prime denotes the derivative with respect to y. 

The horizon is a regular singular point of the differential equation, and demanding 
the regularity of the perturbations on the horizon, an appropriate boundary condition is 
required at the horizon: 

M 1 ) + ( 2 ~ 5=3) M*) = <>• (3.9) 
According to the discussion in the beginning of this section, we take 

6i,o(l) = 2. (3.10) 

From the zeroth law of thermodynamics, all Kaluza-Klein modes of a(x, y) must vanish 
at the horizon. So we have 01^(1) = 0, as is also evident from ( |3.8| ) and the finiteness of 
6i,o(l). 



4 Another type of single linear second-order ODE was discussed in j3fj], based on perturbations of Eu- 
clidean Schwarzschild black holes. However, since the equation has a singular point between the horizon 
and the infinity, it is clearly less tractable compared to our master equation. 
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D = d+1 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


k 


0.876 


1.27 


1.58 


1.85 


2.09 


2.30 


2.50 


2.69 


2.87 


3.03 


3.19 


h 


1.62 


2.20 


2.71 


3.09 


3.29 


3.23 


2.91 


2.29 


1.30 


-0.05 


-1.83 



Table 1: Gregory-Laflamme critical wave number fco for the effective dimension D = 5 ~ 15. ki is 
the next-order correction to the critical wave number. The horizon radius of the uniform solution 
is y = 1. An approximate algebraic equation of the GL critical wave number for D < 9 is given in 
Eq. (l3l4l). 



The perturbation equation has two independent solutions, and at infinity the respec- 
tive solutions are exponentially growing or decaying. The asymptotic flatness implies the 
growing mode is absent at infinity. The asymptotic behavior of the decaying mode is 

ho ~ e~ koy y (d ~ 4)/2 . (3.11) 



Thus, finding the solution of Eq. (|3.6| ) is a one parameter shooting problem, which is easy 
to carry out. The GL critical wave number kg is determined for several dimensions and 
is summarized in Table [l]. 5 The values of ko do of course accord with previous results 
|H, ||, [l(], |3f| . The numerical plots of the solutions are given in Fig. |]. 

Let us take a closer look at the solution near the horizon. We assume the following 
series expansion for b\ o, 



oo 



^/?n(y-i) n . (3.12) 



"1,0 

n=0 

Substituting this into Eq. (|3,6j), we obtain a five-term recurrence equation. The recurrence 
equation can be solved by employing the continued fraction method, which determines the 
critical wavelength. However, the shooting method is more tractable than the fraction 
method, and thus we refrain from presenting the latter. The first few terms are still useful 
to observing the near horizon behavior: 

A = 2, 



ft = 2 ( t0 _ 2+ _M_), 

^(fc + 4d-12) 



p 2 = 8 -2d+^ 



8 + 6k 



(d - 3)2 

/% = •••• (3.13) 

If we truncate the series expansion at some level and match the truncated solution to the 
asymptotic solution at an appropriate radius, we might be able to obtain an approximate 
solution of the critical wavelength. For the truncation at 0{(3^) and matching at y ~ 2, we 



5 Note that we have performed the calculations up to D — 20, although we show only a part of them in 
the Tables. 
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Figure 1: The numerical plots of metric perturbations for D = 6. The event horizon is located 
at y = 1. Each line is normalized by the absolute value of b at the horizon. Only &o,i (and ao,i) is 
the zero mode, so it shows slow decay at the asymptotics, compared to other KK massive modes. 



find a fourth-order algebraic equation 
$ + (d - 3) 

;i38 + (l + d 

d-12 d-12 



4 ,3 {138 + (l + rf)[-49 + 3(l + rf)]}/c 2 4(d-16)(d-3)fcp s 

' j 1 o Jin \ / 



« 0. 
(3.14) 



A real solution of this algebraic equation provides a good approximate solution of the GL 
critical wave for T> < 9 with small errors (a few percent). A more accurate and readily 
available algebraic equation could be obtained by pushing the truncation level up to 0(/?s). 
But the algebraic equation becomes eighth order and it has no analytic solution in general. 

We have seen the existence of the GL critical wave number ko for an arbitrary value of 
a, which does not appear in the master equation (|3.6|). In other words, the smeared black 
branes suffer from the GL instability, irrespective of their charge. As mentioned in Sec. 



2.2, however, the non-extremal smeared black brane cannot be locally thermodynamically 
stable for any value of a. Therefore, the existence of the GL critical wave number does not 
imply the breakdown of the GM conjecture. 

3.2 Second-order perturbation 

The first-order perturbations determine the GL critical wave number ko, and they are 
localized near the horizon due to the rapid decay at the far region. The localized modes 
produce sources for higher-order perturbations that give back-reaction to the background 
geometry. At the second order 0(e 2 ), two types of mode appear: the KK modes, which 
decay exponentially, and the zero modes, which have power-law decay at the asymptotics. 
Since the power-law decay of the zero mode gives finite contributions to the asymptotic 
charges, the zero modes have rather different characters compared with the KK modes. 
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3.2.1 Zero modes Xo,i 

The equations for the zero modes are rather complicated. The ODEs for ao 1 and 60,1 are 

3y 2 - 1 4 

b 0,^ + ^2 _ jT b 0,l ~ y 2 _ ^ (^0,1 - &0,l) = S r 5 ,i(&l,0;feo), 

°o,i + ~T2 — TT a o,l 2~~r( a o,i ~ ft o,i) = S'ao.i^i.o;^), (3.15) 

y{y - 1) ' ^ - 1 

where we abbreviate the source term to 5x 1 which consists of a quadratic expression in 
b±o an d its first and second derivatives. All second derivatives are eliminated from the 
source term using the equation of the first-order perturbations. The critical wave number 
appears in the source term. 

To have decoupled equations, it is convenient to introduce new variables, 

co,i = &o,i - 00,1, d ,i = b ,i + oo,i< (3.16) 

Then the equations are 

„ f + d-3 , 2(d-3) 2 

C 0,l "i 7 c 0,l o7 C 0.! ~~ l! 'co,il & l,0)«oJ) 

yj ' y j 

, f + d-3 j , 2(d-3)(d-l) 
d o,i + ^7 «o,i + ^pj Co - 1 = "Vi (01,0; «o), (3-17) 

where 

Vi( & i,o; fcd) = ~~^~2j i 2 ( d ~ 3 ) 2 ( a i,o - &i,o) 2 + (d - 2)k 2 ,y 2 b 1>0 [2o ll0 + (d - 3)6 1)0 ]} , 

(&i,o5 *b) = -jt {2(d - 3)(d - l)(ai, - 6i,o) 2 + (d - 2)k 2 y 2 b lfi [2a 1)0 + (d - 3)&i, ]} • 

(3.18) 

Note that from the Einstein equations we also have an additional equation which takes 
schematically the form d' 01 = F(&q,1j co,i, b' 1 , c' x ; &o). It is not an independent equation; 
the second-order equation ( |3.17 ) for do,i follows from the first-order equations. 

If we assume the regularity of the perturbed quantities, we obtain the horizon boundary 
conditions by setting y = 1: 

0^(1) = 2{d - 3)00,! - 1 [2(d - 3) + (d - 2)A; 2 ] 6 2 i0 , 

d 0>1 (l) = -2(d - l)c ,i + \ [2(d - 1) + (d - 2)kl)} b\ fi . (3.19) 

Since at the asymptotic region the source terms decay as O(e~ 2koy ), we can neglect the 
right-hand side of ( |3,17j ). Then the leading asymptotic behavior of do,i and co,i are do,i ~ 
l/y( d ~^ and co,i ~ l/y 2 ^ d ~ 3 ' . Interestingly, the power-law decay of co,i is so rapid that it 
does not contribute to the asymptotic charges, due to Eq. fl2.6|) . 

Now let us discuss how we can solve these equations. The shooting parameter of this 
system is co,i(l). The value of do,i(l) at the horizon is not a shooting parameter. It is an 
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arbitrary constant. This is because do,i has a shift symmetry in the differential equation 
and the boundary condition, 



do i — > do i + const. 



(3.20) 



Clearly, we cannot choose the constant in Eq. ( |3.20| ) arbitrarily: If we shift the value of do,i, 
the value of ao.i on the horizon changes, which means that the temperature and entropy of 
the black brane changes. The constant should be fixed by imposing the asymptotic flatness 
on g?o,i, after one obtains a solution. 

The numerical plot of the solutions are presented in Fig. |], in which one can see the 
power-law decay of the solution. If we expand the constant x i n Eq. fl2.6|) as 



x = Xl e 2 + 0(e A ), 

the asymptotic behaviors of ao,i and 6q,i are obtained from do,i as 

Xi 



o.i 



y 



d-3 ' 



(3.21) 



(3.22) 



It is worthwhile to note that from the first law of the thermodynamics the following relation 
holes, 



Xi 



°o,i(l) 



(3.23) 



The results of numerics for various effective dimensions T> are shown in Table ||. 
3.2.2 KK modes X 2 ,o 

The KK modes at this order are 02,0 and 62,0- The function 02,0 can be solved algebraically, 

_ 2f(d - 2){yb 2 ,o)' ai,o [ai,o - (d - l)bi, ] (d - 2)fb lfi 

a2 '° d-3 + (d-l)/ 8 4[ d _ 3 + ( d _i) / ]2 

2(d -3)(d- 2)ai >0 + &i,o {2(d - 2)^y 2 - (d - 3) 2 + / [5 + (d - 4)d]} ) . (3.24) 

Here the source term consists of &i,0;&2,o (and also a^o). A master equation for 62,0 is 
similar to the first-order perturbation: 

62,0 + Ub' 2fi + V(i2)&2,o = (&i,o; ^0), (3.25) 
where the source term #;, 2 is quadratic in 6i ; o and its first derivative. The explicit form is 



i>2,0 



d-3 



a i,o + a i,obi,o< (d - 5) 



2h 2 



4/(d-3)(d-2) s 
[d- 3+ (d -!)/]=» J ^-3 + ^-1)/]^ 



i (d - 3) [3 - d - 2(d - 2)k 2 y 2 } - f {d[13 - 6d + d 2 + 2(d - 2)A;gy 2 ] - 12} 
/ 2 



d- 3 L 



33 - 56d + 36d^ - 10d a + d 4 - 2(d - 2)(d - l) 2 k^y 2 



(3.26) 
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D 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


On 1 ( 1 ) 


-1.63 


-1.30 


-1.20 


-1.15 


-1.13 


-1.11 


-1.10 


-1.09 


-1.08 


-1.07 


-1.01 


60,1 (1) 


-0.65 


0.36 


1.05 


1.65 


2.21 


2.75 


3.29 


3.8 


4.3 


4.9 


5.4 


62.0(1) 


-2.83 


-3.96 


-5.09 


-6.22 


-7.36 


-8.51 


-9.66 


-11 


-12 


-13 


-14 


Mi) 


19.7 


21.9 


20.2 


9.91 


-1.40 


-8.25 


-17.4 


-27 


-25 


-60 


-69 


mi 


1.45 


3.03 


4.69 


7.30 


7.41 


7.95 


7.66 


6.3 


3.6 


-0.67 


-6.8 


si 


2.89 


4.55 


6.25 


7.78 


8.90 


9.27 


8.76 


7.1 


4.0 


-0.74 


-7.4 




-1.04 


-1.08 


-1.11 


-1.09 


-1.01 


-0.85 


-0.62 


-0.31 


0.09 


0.56 


1.1 



Table 2: Summary of numerical results for the static perturbations. The range of effective space- 
time dimensions is T> = D — p = 5 ~ 15. The amplitude of the linear perturbation is set as 
&i,o(l) — 2. The quantities fii and si change their signs above T> = 13, whereas t\ changes its sign 
above D = 12. 



where ai in the source term can be eliminated by Eq. (|3.8|) . 

The asymptotic behavior of these KK modes are &2,o>«2,o ~ exp(— 2fcny), and the 
horizon boundary conditions are 

&i,o (3-27) 

and d2o(l) = 0. Since the critical wave number k§ is given at the first order, this ODE 
is a shooting problem with one shooting parameter, i.e., amplitude at the horizon 62,0(1)- 
Applying the same numerics as for the first order, we can easily obtain the solution. The 
parameter that gives a regular solution is summarized in Table 

3.3 Third-order perturbation 

As we see later, the second-order perturbations are not enough to determine the difference 
of the entropy between a non-uniform black brane and a uniform one of the same mass. 
Knowledge is required of the correction to the critical wave number, ki, which appears at 
this third order. At the third order 0(e 3 ), new KK modes X± 1 and A^o appear indepen- 
dently. The equations of the latter modes do not contain k\ , and so we are not interested 
in these modes. The KK modes X\\ are the higher-order counterparts of Xx t o] the form 
of the homogeneous equation is identical to Eq. ([O]). Then the master equation for b% t i 
takes the form of 

£a&i,l = - 2fe ° fcl 61,0 + S bl l (b Q>1 , bi ;0 , 6 2 ,o; fco). (3.28) 

The function ai,i(y) can be solved algebraically as before. The source term Sb 11 is in- 
dependent of k\, and the algebraic equation for a\ t i does not explicitly contain k\. The 
dependence of a\ t i on the critical wave number k\ is only through b\ t i. The horizon bound- 



62,0(1) + 2 



fi, 



^2,0 



1 (d-2)k$ 

2 d-3 
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ary condition is given by 



2d — 6 — fcg 2k\k bifi 



d-3 1,1 d-3 ' 8(d-3) 

x {-2(d - 3)(8a ,i - 8&o,i + &?,o " 46 2 ,o) + (d - 2)fcg [86 ,i + (d - 2)^ + 46 2j0 ] } , 

(3.29) 

where all functions are evaluated at y = 1. We can solve this second-order equation 
in the same manner used for the first-order perturbations. The main difference is that 
the present case is a two-parameter shooting problem; is not arbitrary. A regular 

solution exists only for an appropriate set of b\ i(l) and fei. We have performed the two- 
parameter shooting, and the results for various dimensions are summarized in Table [l] and 
|2] (see also Fig. p]). We should notice that the sign of k\ changes above T> = 13: 

h > (D < 13), 

ki < (D > 13). (3.30) 
This would be a precursor to the existence of a critical dimension. 



4. Thermodynamics 

We are now ready to calculate the thermodynamical quantities of the perturbative non- 
uniform solutions. From Eq. (|2.5| ), we obtain the mass M, temperature T, entropy S 
and charge Q for a non-uniform smeared black brane. The differences for these quantities 
between the non-uniform and the uniform critical solutions are 

*JL = ^ + o(^) 

M (d-2)/(d-3) +sinh 2 a V ; ' 

ST ' l = _I a0)l( l )e 2 + O(e 4 )) 



T VAh 2 
^ = y /T h -l = l a0 , 1 (l)e 2 + O(e% 

f-a 

These physical quantities are expressed in terms of the perturbative quantities, which are 
calculated without fixing the asymptotic periodicity of the circle. It is easy to compare the 
obtained physical quantities with others if they are given in a frame where the asymptotic 
length of the circle is fixed. Without transforming physical quantities into such a frame, 
we introduce variables that are invariant under rigid rescalings of the entire solution. Once 
physical quantities are expressed in terms of the invariant variables, we do not need to care 
about the variation of the asymptotic periodicity. 
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The invariant quantities can be obtained by multiplying K by suitable powers: 

8n_8M 5K 2 4 

7 = ^- + (d-3)— = W£ +0(e ), 

5t 5T 5K 2 4 

£ = f +(d-2)f =. ie > + 0(A 

| = f + (d-a)^ = « £ 2 + o( £ 4 ), (4.2) 

where the leading-order coefficients are given by 



/M (d-2)/(d-3) + sinh 2 a 1 j fc ' 



1 . . k x 
n = -o a o,i(l - — , 

2 fe 



*1 = ~a ,i(l) + (d-2)-i, 
2 fc 

gi = (d-3)£. (4.3) 

Here we see the reason why we have integrated the perturbations up to order 0(e 3 ). The 
correction to the wave number k\ is necessary to obtaining invariant quantities. 
Note that non-extremal solutions obey the first law, 

dM(R ,a) = T(R ,a)dS(R ,a) + <5>(a)dQ(R ,a), (4.4) 

provided the asymptotic circle length is fixed. We can confirm this directly from Eq. ( |2.5| ) 
by utilizing the first law for the uncharged system, which gives 

1 dA h d-3 dx 

X--77T- ( 4 - 5 ) 



2A h dR R dR 
4.1 Vacuum black branes 

To begin, let us first consider the thermodynamics of neutral black branes for vacuum 
spacetime. 6 What we are interested in is the difference between the entropy of a uniform 
brane and that of a non-uniform one of the same mass. The difference is evaluated as 
follows: 

SnU — Sy _ 2 4 . 6 , 



Sv 



aie 2 + a 2 e 4 + 0(e t 



d-3 d-2 ( 1 \ , . 

ai = ^ ~ d^2 S1 ' a2 = -2{d^3){ Tl + d^3 f11 )^ (46) 

where we have used the Smarr formula TS/M = (d — 2 — n)/(d — 1) and the tension of 
uniform strings, n = \/{d — 2). This equation can be easily derived by starting from the 



6 Since there is no charge in this system, the metric ( |2,3| ) is identical to that of the (neutral vacuum) 
black (p + l)-brane. 
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frame with SK = and re-expressing the results in terms of the invariant quantities ( |4.2j ) . 
See Appendix [B| for more details of the calculation. 

As discussed above, the relation ( |3.23| ) holds for \i an< l 01,0(1) from the first law, 
and then o\ vanishes. As a result, the entropy difference arises at 0(e 4 )- The correction 
02 becomes 02 = cto.iMi/ '4. Since the numerical value of ao,i is negative (Table |2|), the 
sign of o"2 is determined by /ij. Our numerical results for D = 4 ~ 15 are summarized in 
Table |2[ and they are consistent with previously obtained results [jl], |3(| ; 02 is negative 
for T> < 13, hence the non-uniform phase is not entropically favorable and it implies a 
first-order phase transition. For T> > 13, 02 becomes positive, implying a higher-order 
phase transition as observed in Ref. 36]. These are summarized as 

02 < (D < 13), 

02 > (D > 13). (4.7) 

The vanishing of o\ can be used as a check of numerics. In fact, our numerical re- 
sults give vanishingly small a\ for relatively low dimensions (D < 9). For much higher 
dimensions, however, it is hard to read off xi numerically because of the rapid decay of 
the zero modes. Thus, we have calculated xi by using the relation ( |3.23| ), which allows 
very precise determination of physical quantities because ao,i is a local quantity. Moreover, 
the results obtained here are very simple and concise compared to those obtained by the 
general perturbation theory ||; the invariant quantities (4.2) under rigid rescalings are 
given only by 3 variables, i.e., ao 1, &o> an d k±, whereas for the general perturbation theory 
many quantities enter into the estimation of the scheme independent quantities. 

Next, we consider a canonical ensemble, i.e., fixed temperature. Helmholtz free energy 
F = M — TS for the neutral solutions is given by 

^ = sf^V 3 [l-(^-3)x]. (4.8) 



n 2vr r 

The free energy of the uniform strings is F\j = 9[(d — 3)/A7rT] d ~ s . The difference of the 
free energy between the uniform and the non-uniform solutions for the same temperature 
is evaluated as 

' Nl h P2e 4 + 0(e 6 ), p 2 = ~{d- 3)(d - 2) ( n + -^—^1) n. (4.9) 



Fu 2 V /v ' \ d- 3' 

See Appendix [B] for the derivation of Eq. ( |4.9| ). The structure of the coefficients P2 is quite 
similar to that in the microcanonical ensemble ([O]), except the overall coefficient, which is 
n in the present case. It is interesting that the invariant variable n changes its sign at the 
different dimension D from that for n\. t\ is negative for D < 12 and positive for T> > 12 
(Table |2[). Consequently, the free energy of the non-uniform phase becomes favorable at 
D > 12; 

p 2 > (D < 12), 

P 2 < {D > 12). (4.10) 

It is remarkable that T> = 13, which is the critical dimension in the microcanonical ensem- 
ble, is not special in the canonical ensemble. 
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4.2 Weakly charged smeared black branes 

General thermodynamic properties of the charged (dilatonic) black branes are rather com- 
plicated compared with the vacuum cases. Besides, the general formula that is given in 
terms of the invariant quantities is untractable, so we focus attention on two limiting cases: 
a weakly charged case and a near-extremal case. The latter will be discussed in the next 
subsection. 

Firstly, we consider weakly charged smeared black branes in a microcanonical ensemble 
(fixed mass and charge). Expanding the entropy difference for small charge Q <C M, we 
find 

2 

5o 2 + 0(Q A ) 



Snu — S\j 



Sv 



(72 + 



Q 



M- 



Sao 



(d-2) 



Ti + qi 



2(d-4) 
d-3 



"Mi 



qi + 



d 



d-3 



Ml - n Mi 



(4.11) 



2(d-3) 2 

where p± and n are evaluated at a = 0, and we note that Q/M is invariant under rigid 
rescalings. In Appendix ^, we will show how Eq. ( |4.11| ) is derived. As in the vacuum case, 
the leading entropy difference appears at 0(e 4 ), and no correction appears at a lower order. 
This is not due to an accidental cancellation under the approximation of weak charge, but 
we can show after a tedious calculation that leading corrections are always 0(e 4 ) for general 
charge Q. The correction term 5<J2 is positive with numerical values of 0(1) for D < 15, 
whereas 5a2 becomes negative for a very large value of T). (We have confirmed this up to 
T> = 20.) In any case, the new correction is sufficiently small compared to 02, as long as the 
assumption of weak change is valid, and hence the charge and dilaton do not change the 
thermodynamic phase structure. One might think that the phase structure is drastically 
affected by the bulk fields with an arbitrary charge, because the correction 802 counteracts 
the leading term 02 in the above case. As we will see below, however, the same result of 
Q4.7| ) holds even for the near-extremal limit in a microcanonical ensemble. 

Next, we turn to a canonical ensemble for weakly charged branes. The difference of 
the free energy for the same temperature and charge is 

+ ( Q \\ ~ 

P2 + T7 <V2 



Fnu — Fu 



Fu 



M 



SP2 



(d-2)' 



, )((] . J(d - 3)( ai + n)ri - [qi + 2{d - 3)n] (4.12) 

See Appendix || for the derivation of Eq. fl4.12| ) . The dimensional dependence of 5p2 is 
depicted in Fig. ||. From the figure we see that 5p2 is negative for D < 13, hence, the 
correction counteracts P2 ■ The correction in this weakly charged case is too small to change 
the phase structure; interestingly, the effect becomes evident near the extremality, as we 
will see in the next section. 



4.3 Near-extremal smeared black branes 

Let us finally consider the thermodynamic behavior near the extremal state. A near- 
extremal limit can be realized by taking a to infinity. We define energy above extremality 
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Figure 2: Dependence of the coefficients on the effective spacetime dimension D. p^ and 8p2 are 
divided by D and D 2 , respectively. The sign change of 02 is at D > 13, whereas the sign change of 
p 2 is at D > 12. 



M, "reduced temperature" T and "reduced entropy" S as follows: 



7= lim JQT 

a— »oo 

s 



(rf-3) 3 S (LRq\ 
AA h I?B* V 2vr ) 



d-3 



1/2 



8 = 

M E 
M e i 



lim 



l67r 2 A h S ( LRq 
(d-3) 



d-l 



lim(M-Q) = g(^o 
= lim (M ei - Q) = - 



2tt 

d-3 



1/2 



d-l 



(d - 3) X 



LR \ d ~ 3 d-3 



(4.13) 



2vr / 2 

Note that M — M G i = M\ a= Q and T8 = TS"| a= o hold. Using Eq. ( |4.5| ), we can rewrite the 
first law in terms of the reduced quantities near the extremal state: 



dM(R ) = 7(R )d§(R ). 



(4.14) 



In addition, we notice that the Smarr formula also takes the same form as that for the 
uncharged case. 

The invariant fractional changes of the reduced quantities are given by 

We 2 + 0(e 4 ), 

:S - ie 2 + 0(6 4 ), (4.15) 



5fi 


8M 




3>£ 






6f _ 
f 


T + 






6s _ 
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<5§ 

Y + 
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V 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 




-0.53 


-0.91 


-1.32 


-1.71 


-2.00 


-2.10 


-2.00 


-1.6 


-0.9 


0.3 


1.9 


P2 


0.18 




-15.7 


-10.4 


-9.19 


-8.40 


-7.43 


-6.0 


-3.8 


-0.8 


3.3 



Table 3: Numerical results of invariant quantities for near-extremal smeared black branes. 



where the leading-order coefficients are given by 

2(^-3) fci 

" 1 = — d=r xi + (d - 3) v 
i i , &i 

n = --a ,i(l) + -( a -5 )^' 

S -i = ^a 0il (l) + i(d-l)|. (4.16) 

We consider first the thermodynamics in a microcanonical ensemble (fixed energy M). 
The entropy difference between the uniform and non-uniform solutions for the same energy 
M. is computed with the first law ( 4.14] ) as 



§NU — §U = 4 , ^,Jn = d—l 



cx 2 e 4 + 0(e b ), a 2 



4(d - 3) 



d-5 
2(^3) W " Tl 



Mi. (4.17) 



The numerical coefficients of /xi and t\ differ from those of ( f4.6| ), due to the "corrections" 
M e i to the mass in the Smarr formula (gOD , i.e., 18/ (M - M e i) = 2(d - 3)/(d - 1) for the 
uniform phase (x = 0). The entropy difference again arises at 0(e 4 ) as the neutral case. 
As we can see from Table ||, the entropy of the non-uniform phase is less than the entropy 
of the uniform phase for a given energy M for D < 13. In other words, the non-uniform 
phase is disfavored thermodynamically over the uniform phase. For D > 13, on the other 
hand, the non-uniform phase is favorable over the uniform phase: 

a 2 <0 (D < 13), 

ct 2 > (V > 13). (4.18) 

The critical dimension at which the sign of a 2 changes is same as the one for the uncharged 
case. 

We now turn to a canonical ensemble (fixed temperature 7). We define the "free 
energy" of the solution as 

7 = M - 7$. (4.19) 
The free energy of the uniform solution is 

7u(Ro) = -^M(i? ), (4.20) 

which vanishes for d = 5 and is negative (positive) for d > 5 {d < 5). The difference of free 
energy between the uniform and non-uniform solutions for the same temperature is given 
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by 



^nu -3tj - 4 ,/->/■ 6\ - d-l (I _ d - 3 _ , 

/>2£ + 0(e ), P2 = 3 — r -Mi - - — -n ti, (4.21 



where D / 6 [d ^ 5). The exceptional case of D = 6 has to be treated separately because 
the free energy of the uniform solution vanishes. This case will be discussed elsewhere 
(see Ref. p7|] ). Table || summarizes our results. We see that p2 is positive for D = 5 and 
D > 14. This change of sign at the lower dimensions is in contrast to the previously studied 
cases in which the change of sign occurs only at relatively higher dimensions. However, we 
should recall that the free energy 3~u of the uniform solution changes its sign depending on 
the effective dimension. Consequently, the non-uniform phase becomes thermodynamically 
favorable over the uniform phase only at T> > 14: 

Ju<^nu CD<14), 

Ju>^nu (D>14). (4.22) 



This result differs from the vacuum case ( 4. ldf ) . We have observed in the weakly charged 



case that the correction due to the charge has a tendency to change the phase structure. 
The large correction near the extremality shifts the critical dimension upward, although 
the correction cannot change the phase structure significantly. 

5. Conclusion 

We have investigated the uniform and non-uniform black branes smeared on a transverse 
circle by constructing perturbative solutions explicitly. We have made use of the HO metric 
ansatz in our perturbation analysis. In the original perturbation scheme Q based on the 
general metric ansatz, analysis requires a good deal of subtle maneuvering. Making use 
of the HO metric has made the analysis enormously simple and clear in many respects. 
First of all, the EOMs of the charged dilatonic black branes are the same as those of the 
vacuum black branes on a circle. Consequently, we do not need to perturb the dilaton 
and form fields at any order, and the task is reduced to solving only one master equation 
at each order. Owing to such simplicity, the asymptotic behavior of the perturbation and 
the (analytically approximate) GL critical wave number are easily obtained. Besides, all 
physical quantities are characterized by a single parameter ao,i- 

The existence of a static mode means an onset of unstable modes, hence, the smeared 
black branes suffer from the GL instability irrespective of their charge. However, this 
does not mean that they provide a counterexample of the GM conjecture since there is 
no parameter region in which local thermodynamic stability is realized, except at the 
extremality. 

Having performed the perturbation up to the third order, we obtain the phase structure 
around the critical uniform solution. We have begun with studying vacuum black branes. 
In the micro canonical ensemble, the thermodynamic system has the critical dimension 
D* = D — p = 13, at which the nature of phase transition between the uniform and non- 



uniform branes changes [36]. On the other hand, we found that the critical dimension 
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changes if we change the ensemble. In the canonical ensemble, the critical dimension is 
D* = 12. Furthermore, if the system is near the extremal state, the critical dimension in the 
canonical ensemble becomes 25* = 14, whereas the critical dimension in the microcanonical 
ensemble does not change. We think that a specific value for a critical dimension is not 
universal. It depends on matters in the bulk and type of ensemble. 

While the third-order perturbations have been performed to obtain meaningful leading- 
order corrections, the next-order corrections are also interesting to study. To obtain such 
corrections, we need to perform the perturbations up to the fifth order, and it would 
appear possible to obtain them. We think that the perturbative results in this paper 
give a good theoretical basis for understanding and constructing fully nonlinear solutions 
of the HO ansatz. Fully nonlinear solutions will allow us to explore the thermodynamic 
phase structure for an arbitrary charge, which might have significant consequences. We 
will discuss the fully nonlinear (numerical) solutions in a forthcoming paper. It is clearly 
interesting if a critical dimension appears at much lower dimensions, because there is no 
sensible quantum theory of gravity for D > 10, 11. Such a situation might be possible if 
we consider other types of bulk fields or configurations. 
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A. Vacuum Einstein equations 



In this appendix we give explicit representation of the Einstein equations. From the ansatz 
( |2.3| ), the metric of the (d + l)-dimensional vacuum black branes on a circle is 



ds 



d+l 



e d tf + e «-id-2)b dv 2 +ebR 2 d tf 
f 



(A.l) 



where we have taken A(R,v) = e a and B(R,v) = e b . The nontrivial components of the 
vacuum Einstein equations ft*- = are given by 



rnv 
Ji,. 



cr>t) 



— (f) (/>/ 2)//' a" r- 1 



+ 



(d- 2)6 + a 

I. | 2 (d - 3 + /) [(d- 2)6' - a'] - e^ b R(d - 2)(d - 1)6 2 } 



+(d-2)dbe (d - 2)b + (d-2)f 



2 , (d-l)b> 

a H 

R 2 



b'- a - 
R 



er a /27rV 



2(d-3) b 
Ri [ ' 



- 3 + f)b' + e^ b (d - 2)6 2 + fb" + e^ 2 H 



R 
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R 



T"(t) (^{ d " 3 + /[ d - 1 + (' i - 2 )^ / ]} 



+(d-2)/ 



- - (d - 1)6' 



b-2f{d-2)b'\, (A.2 



where the prime and dot stand for the derivatives with respect to R and v, respectively. 
Since the angular component Jig does not contain any derivatives of a(R,v), a(R,v) can 
be given in terms of b(R, v) and its derivatives. 



B. Entropy and free energy difference 

In this appendix, we derive the formulae of the entropy difference and free energy difference, 
i.e., Eqs. ([O]), flOJ), ( pL"TTD and ( |4T2|) . 

First, we expand the thermodynamical quantities by e as follows: 



M = Y^M p e 2 ^ T = Y,T p e 2 ? 1 S = ^S P e 2p , <? = £<? P e 2p - (B.l) 

p=0 p=0 p=Q p=0 

The entropy difference between the uniform and non-uniform branes of the same mass and 
charge is expressed as 



5*nu — Su Sp=o S p e 2p 



1 



__ 5i 2 S2 4 _ d\ Agp 2 _ Affp / A5 \ , 6 ^ 

Dq Do Do Do Do \ Dq ' 



(B.2) 



where ASo is the entropy change of a uniform brane due to the change of the mass and 
charge, which we denote by AMo and AQq, respectively. 

Let us focus on the quantities Si /So and S2/S0 in Eq. (|B.2| ). From the first law with 
a fixed scale of the circle (5L = 0), we have 



Mi = TqSx + Qi tanh a, M2 = T0S2 + ~^TiSi + Q2 tanh a. 



(B.3) 



From Eq. fl2.5|), the mass of a uniform brane is written in two ways: 



d-2 
d-3 



+ sinh a) To So and Mi 



1 + 



d-2 



(d — 3) sinh a 



Qotanha. (B.4) 



From these equations flB,3|) and (|B.4| ), we can write the entropy difference, Si/ So and 
S2/S0, in terms of Mi/M , M 2 /M , Q1/Q0, Q2/Q0 and T\/Tq. 

Next, we focus on the quantity ASo/So in Eq. ( p3 . 2|) . To express the entropy of uniform 
brane in terms of mass and charge, we solve Eq. ( |2.5| ) for a and Rq with \ = 0, 



sinh 2 a 

LRp 
2vr 



2(1 - Qg /M 2) 
S(d-3) (d-2 



-1 + 2 



M 



d-3 



+ sinh 2 a 



d-2 /Qo 
d-3 \M 

l/(d-3) 



+ Wi + 



4(d -2) (Q, 



(d-3) 2 \Mt 



(B.5) 
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Substituting above two relations into the entropy of the uniform brane (A^ = 1) in Eq. (|2.5| ) , 
one obtains in the leading order of Qq/Mq (^C 1) 

.S,(.U () .<5o)«47rS .. il ' 



r M 


(d-2)/(d-3) 


L(d-2)S_ 





d- 3 



M 



(B.6) 



From this expression, we can compute A So = S(Mq + AMq, Sq + ASo) — S(Mq, Sq). Taking 
the mass and charge differences as AMo ~ M\e 2 + M2e A and AQo ~ Qi^ 2 + Q2^ 4 , we obtain 
ASo/So in terms of Mi/Mq, Qi/Qo and so on. 

Finally, replacing such gauge dependent quantities with gauge independent quantities 
such as Hi with 5K = 0, we obtain the expressions (4.6) and ( 4.11| ). 

In a similar way, we can derive the formula of the free energy difference, Eqs. ( f4.9| ) and 
(4.12). The free energy difference between uniform and non-uniform branes of the same 
temperature and charge is given by 
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To F F 

where AFo means the free energy change of a uniform brane due to the change of temper- 
ature and charge. From the first law with a fixed scale of the circle, we have 



1 



-5 T 2 + 2 ^i^i + ^2 tanh a. 



From Eq. (| 
F 



Fi = —SoTi + Qi tanh a, F2 
), we can express the free energy of a uniform brane in two ways: 
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(B.8) 
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From Eqs. ( |B.8| ) and ( |B.g| ), we can write the free energy difference, Fi/Fq and F2/F0, in 
terms of Ti/Tq, T 2 /T , Qi/Qo, Q2/Q0 and Sb/Si- 

To calculate AT0/T0, we have to express the free energy of uniform branes as a function 
of temperature and charge. From Eq. ( |2.5| ), we have 
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+ 0(Qg). 



With these relations, the free energy of a uniform brane is given by 



T (T ,Qo) = S 
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(B.10) 



(B.ll) 



From this, one can compute the free energy difference of a uniform brane as ATo = To(To + 
ATo,Qo + AQo) — To(To,Qo)- Taking the change of temperature and charge as ATo « 
Tie 2 + T 2 e 4 , AQq « Q\e 2 + Q2£ 4 , we can write AFq/Fq in terms of Si/Sq, Ti/Tq and so 
on. Therefore, replacing these with gauge independent quantities, we obtain Eqs. (| 
and flOp . 
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